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Abstract 

The zero momentum limit of thermal self-energies calculated in 
perturbation theory depends on the order in which the time and the 
space components of the momentum are taken to zero. We show that 
this is an artifact of the perturbative calculation, and in fact the limit 
is well-defined when higher orders of the perturbation expansion are 
properly taken into account. 

The existing calculations of thermal self-energy functions using the for- 
malism of Thermal Field Theory yield results that are not defined if the 
external momentum 4-vector is zero. The classic example is the photon 
self-energy vr^,^ in an electron gas. It is well known [|l| that the result of 
the one- loop calculation of T^^uiv^ iP) for a photon with external momentum 
pfj. _ (jP^p^ is such that 

lim ^^1.(0, p) ^ lim ^^^(/,0) , (1) 
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so that the hmit in which all components go to zero is not defined. To be 
more specific, the above-mentioned problem occurs for the real part of the 
self-energy, while the imaginary part is well defined. 

Various attempts have been made to resolve this puzzle 0, |^, H, |^, which 
involve either introducing new and ad-hoc Feynman rules for thermal field 
theories, or putting restrictions on the general rules. Along another line of 
attempt, it has recently been pointed out by Arnold, Vokos, Bedaque and 
Das (AVBD) |^ that the problem mentioned above occurs only if the self 
energy diagram contains two propagators of the same mass. If the masses of 
the particles in the loop are different, the problem does not exist. In fact, 
the calculations of the neutrino self-energy in a gas of electrons and nucleons, 
which were carried out even before the work of Ref . ^ , show this [0, |§, . It 
has been speculated that this property may be utilized to introduce a mass- 
splitting regularization for thermal diagrams |^ in cases where problems are 
known to occur. 

This problem, as well the attempts to resolve it, are based on the re- 
sults of one-loop perturbative calculations. It is natural to ask whether the 
singular behaviour of the self-energy function at zero momentum might be a 
consequence of the approximations and idealizations that are implicitly made 
in the perturbative calculations. In this article we show that this is precisely 
the case, and that the calculation of the self-energy beyond one-loop order 
yields a result that is defined at zero momentum. This result is obtained by 
calculating the self-energies using the full propagators for the particles that 
appear in the internal lines of the loop diagrams. Some of these propagators 
have an absorptive part which, as observed by AVBD 0, is well-defined at 
zero momentum even if they are calculated to one-loop and the particles in 
the internal lines of the loop have the same mass. As we will see, this in turn 
governs the zero-momentum limit of those self-energy diagrams in which the 
internal lines have the same mass. We exemplify these assertions by calcu- 
lating the photon self-energy in scalar QED, but similar considerations apply 
to QED proper as well. 

Before proceeding, we recapitulate some results of the canonical approach 
to the thermal propagators, which we will be using throughout [jlO[- In 



this approach, one has to use anti-time-ordered propagators in addition to 
the time-ordered ones, as well as propagators with no time-ordering. These 
propagators can be arranged in the form of a 2 x 2 matrix. For example, for 
any bosonic field where A denotes any Lorentz index carried by the field 
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(none for a scalar field, one for the photon), we can write 

i'D^.^ix-y) ^ lT^\x)^^*{y)) , (2) 



iT>^i{x-y) ^ (T$^(x)$^*(y)) , (3) 
i-Dtiix-y) ^ {^^*{y)^\x)) , (4) 
i-Di^x-y) ^ l^\x)^^*{y)) , (5) 



where the time-ordering and the anti-time-ordering operators T and T are 
defined as 

T$^(x)$^*(?/) = Q{x,-y,)^''{x)^^\y) + e{y,-x,)^^*{y)^''{xm 
r$^(x)$^*(|/) = e(yo-Xo)$^(x)$^*(y) + e(a:o-yo)*''*(|/)$^(^)(7) 

© being the step function. We can now make the momentum space expansion 
of the field as 

^^^"^^ = / ^ H^^^^P^ ^)^"^''' + ^a(p)^^(p, A)e^^-] , (8) 

where u"^ and represent different plane wave solutions arranged by the 
index A, and ax{p) and hx{p) are the annihilation operators for particles and 
antiparticles, respectively. For a self-adjoint field like the photon, ax{p) = 
b\{p). The properties of the thermal bath come in from the expectation 
values 

{ax{p)al,{p')) = {2TTf2ESip-p')6xx'[l + fBip,a)], (9) 
{bx{p)b*y{p')) = {27rf2ES{p-p')Sxx'[l + fB{p,-a)], (10) 

with 

fB{P,(^) = ^fSp.u-a_l ^ (11) 

where a plays the role of a chemical potential. We have introduced the 
velocity 4- vector u'^ of the heat bath, which has components (1,0) in its own 
rest frame. 
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For scalar fields, the procedure described above gives []T0[ the 2x2 free- 
field propagator as 



where 



pZ _ _|_ 



and 



^ 'l + r]B{p,a) \ -^{p-u)fB{-p,-a) l + Vsip^a) /' 



with e{x) = Q{x) — Q{—x), and 

?7b(p, a) = Q{p ■ u)fBip, a) + Q{-p ■ u)fB{,-p, -«) . (15) 
It is straightforward to check that this gives, for example, 

Aii(p) = — ^^—7: - 2Tii5{p^ - m^)r]B{p, a) , (16) 



which is the propagator given by Dolan and Jackiw [|TT[|. The explicit forms 



of the other components are also given in the literature [12 



The one-loop diagrams for the photon self-energy in a background of 
particles are depicted in Fig. |l|. Diagram |l]b produces a term that is 
proportional to the total electric charge of the system. However, since that 
term is a constant, independent of the photon momentum, it will not be 
relevant for our discussion and we will not consider it any further. 

The result of calculating Diagram ^ using the free-field propagator given 
above for the (j) field is 

-» f d^p 
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where the vacuum contribution has been omitted, as we will always do hence- 
forth whenever we write explicit expressions for the self-energies. The above 
formula reveals the problem to which we alluded in Eq. (0) . As we mentioned 
earlier, this problem vanishes if the diagram is evaluated employing the full 
propagator of the field instead of the free-field propagator. The full scalar 
propagator, which we denote by A'(p), can be written just like in Eq. (0) 
but with Aq replaced by 

^0 = \—Fr^ (18) 

where Ho is the self-energy function for the (p field. Thus, 

A'ip)=Un(^^^ -A-)"^' ^^^^ 

and in particular, 

A[,{p)=A', + [A',-A'*]Mp,<^)- (20) 

The important consequence of replacing Aq by Aq is that the (5-function 
present in Eq. ([I6| ) is now smeared if the absorptive part of Ho is non-zero. 
For this reason it is the easy to see that, if Diagram |I]a is calculated with 
the propagator A'^^^ for the (j) field instead of the free particle propagator, 
then the problem at zero momentum does not arise in evaluating [77^,^ (A;)] n. 
Notice that the dispersive part of IIo plays no role in this argument. It is not 
difficult to see that retaining only the dispersive part of IIo and neglecting its 
absorptive part does not remove the singularity of [tt^j^Jh at zero momentum. 

The next step is to calculate Uq and show that in general it has an absorp- 
tive part. To this end, we recall that the inverse of the full scalar propagator 
is given by 

A'-\p) = p"^ - TTi^ ~ U , (21) 

where 11 is a 2 x 2 matrix whose components must be calculated using the 
Feynman rules of the theory. Comparing Eqs. (^) and (PT|), the following 
relations are obtained, 

Uu = no + (no-n*)77B(p,«) (22) 
= -n* + (no-n*)r7B(p,«) (23) 

= -iUo-ni)eip-u)fB{p,a) (24) 
Hai = -(no-n*)e(-p-i/)/B(-p,-a). (25) 
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from which it is easily seen that 



Re Hob) = Renn(p) (26) 
ImUoip) = ^. (27) 

Therefore, to determine Ho we must calculate XIn and H12, which can be 
done by evaluating the diagrams in Fig. ||. Since the dispersive part Ho is not 
relevant for resolving the zero momentum problem of the photon self-energy, 
we will not calculate it here. However, for the consistency of our scheme 
it is important to stress that since the internal lines in Fig. |^ correspond 
to particles of different mass, the function Rellii (and hence Re Ho) does 
not suffer from the zero momentum problem according to the observation of 
AVBD §. 

We now turn the attention to the calculation of the absorptive part of IIo. 
The simplest way to proceed is to calculate and then use Eq. (^). As 
in the case of Diagram |l|b for the photon self-energy, the diagram in Fig. ^b 
is irrelevant for our purpose. In Fig. |^a, the scalar propagator that enters for 
is obtained from Eq. (|T2|) as 

Ai2(p') = -27n6ip'^ - m^)Mp\ a)e{p' ■ u) . (28) 

In addition to this, we also need the thermal photon propagator, whose 
form depends on the gauge that is chosen. Here we will use the Coulomb 
gauge which, for several reasons, has been advocated as a convenient one 
for carrying out finite temperature calculations in QED |14[. The only 



component that is needed for the calculation at hand is the 21 component 
which, borrowing from Refs. [ITl, is given by 



£)^r(/c) = 2m6{k^)fBi-k, 0)e{k ■ u){-Sn , (29) 
where ^ 

SuU = QUU ^ T^kuky :^{Uaky + knlly) . (^O) 



Here uj and n are defined by 



u = k-u, K = V uj"^ — k"^ , (31) 
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and they represent the energy and momentum of the photon in the rest frame 
of the heat bath. A particularly useful relation is 



J2 e^{k,X)ty{k,\) 



A=l,2 



(32) 



where the polarization vectors e^{k,X) = (0,e(A;,A)) are such that, in the 
rest frame of the heat bath, 

k-e{k,l) = k-e{k,2) = 0. (33) 

The application of the Feynman rules to the diagram of Fig. ^ gives 



~^UM = {-ie){ie) J ^^^Di^^^{k){p + p'^A^ip + p'), , (34) 
where we have defined 



k = p' — p . 



(35) 



Substituting the photon and scalar propagators into this expression and using 
Eq. (^71) we obtain 



Imno(p) = --{27i)\{p ■ u) 



(M 



-Mk'')5{p'^ - m')t{p' ■ u)e{k ■ u) 



X {-snip + P'),iP + p'UfBik., 0) - fsip', a)) 
where the identity 

/B(-A;,0)/ij(p',a) = fBip,a)[fBip,a) - /i?(fc,0)] 
has been used. Eq. ( |36D can be written in the form 

d^p' d^k 



(36) 



(37) 



Imno(p) 



-e{p ■ u) 



{2nf2E' {2-Kf2uj 
X { 5^^){p +k-p'){p + p')^,{p + pO.L/, - U] 

+5(4)(p ^k-p'){p + p'),{p + p'Ul + f, + u] 
+^(4)(p + k+p'){p- p'),{p - p\[l + + / ] } , (38) 
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where, for the sake of brevity, we have called the photon and (p particle 
distributions by 

A = /i?(^,o), 

U = fB{p',a), (39) 

while is the antiparticle density distribution, which is obtained from 
by changing the sign of a. In addition. 



p'f' = (E', p') , E' = + m2 , (40) 

k^' = {uj,k), u=\k\. (41) 

A more physically intuitive representation of this formula can be obtained 
by using Eq. ( |52D for Sfj_^. Then, introducing the amplitudes 

Ma = {-ze){p + p're,{k,X), (42) 
Mb = {-ie){p-p're,{k,X), (43) 

we get 

lmUo = -\p-u\T{p) , (44) 

where we have defined 

~ 2p-uJ (27r)32a; {27ry2E' ^ 



X { 6^'\p + k - p')[f,il + U) - Uil + /,)] \Ma? 

A=l,2 

+5^'^{p-k+p')[f^{i+f,)-f,{i+j^)] y: \mb\' 

A=l,2 

+5(4) _ k - p')[(l + + /^) - f^f^] Y \Ma\' 

A=l,2 

+S^'\p + k+p')[fJ^ - (1 + mi +7^)] E \Mb\' } .(45) 



A=l,2 



The formula for T given in Eq. ( ^51) is immediately recognized as the total 
rate for a (f) particle of energy p^ and momentum p (as seen from the rest 
frame of the medium) with integrations over the phase space weighted by 
the statistical factors appropriate for each process |T5[ . Ma is the amplitude 
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for 70 — i> (/) or the decay </> — > 70, while Mb is the amphtude for 00 — >■ 7 or 
700 0. The amphtudes for the inverse reactions are given by the complex 
conjugates of Ma and Mb- For certain specific values of and p some of 
these processes will be kinematically forbidden, but in general F is non-zero. 

In conclusion, the one-loop photon self-energy calculated with the full 
propagator given in Eq. ( pOD instead of the free propagator of Eq. (0), is 
defined at zero momentum. In particular, the absorptive part of the self- 
energy, which physically is related to the damping rate of the particle, cannot 
be neglected if the photon self-energy is evaluated at zero momentum. Then, 
the physical picture that emerges is the following. The traditional formulas 
that are given for 

lim7r^,(0,p) (46) 

IpI^o 

and 

lim7r^,(/,0), (47) 

which are related to well known physical quantities such as the plasma fre- 
quency and Debye radius, are valid in the limiting cases 

/ = ; F < IpI < m 

p = ; F<j9°<m. (48) 

Since the two limits correspond to two different physical situations the results 
are different. Traditionally F is omitted in the above conditions, but then 
it must be kept in mind that the formulas cannot be taken literally all the 
way to zero momentum. The same conclusion can be reached for fermionic 
QED and other field theories, which will be discussed in detail in a future 
publication. 
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Figure captions 

Photon self-energy diagrams. 
Self-energy diagrams for the scalar. 
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